Analytic solutions for nonlinear waves in coupled reacting systems 
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We analyze a system of reacting elements harmonically coupled to nearest neighbors in the con- 
tinuum limit. An analytic solution is found for traveling waves. The procedure is used to find 
oscillatory as well as solitary waves. A comparison is made between exact solutions and solutions of 
the piecewise linearized system, showing how the linearization affects the amplitude and frequency 
of the solutions. 
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I. INTRODUCTION 

Consider a lattice consisting of a reaction system at 
each site n, described by a site potential V(u n ) with u 
a one dimensional field. Neighboring sites in the lattice 
interact through a potential <j){u) . The equation of mo- 
tion: 



dt 2 



-<t>'{u n 



y) + 4>'{un+i-u n )-V'{u n ). (1) 



describes, for instance, the time evolution of the ampli- 
tude of an oscillator, the population of some ecological 
species, or the density of a chemical component of the 
system. If <fi is quadratic, the continuous form of Eq. (|l|) 
is a wave equation augmented by an additive nonlinear 
term: 



d 2 u 2 d 2 u , 



(2) 



where v is the speed of linear waves in the absence of the 
nonlinear potential V. The latter is determined by the 
details of the system under consideration. A choice of 
wide applicability is the logistic reaction term. Equation 
(@) is then: 



d 2 u 



2 9 2 u , 



(3) 



whose interest resides in its relevance to chemical and 
population dynamics, where the reaction term models 
an autocatalytical "malthusian" growth of u, with a sat- 
uration due to intraspecies competition. We have taken 
rescaled variables in Eq. (^) so that the carrying capacity 
of the medium is 1, in order to simplify the notation. 



Equation (||) is very closely related to the so-called 
Fisher's equation [jjj for reaction-diffusive systems, in 
which the left hand side of Eq. (^) is replaced by the first- 
order time derivative du/dt. The two equations share 
the same nonlinear reaction term but, respectively, rep- 
resent opposite limits of completely coherent (wave-like) 
and completely incoherent (diffusive) transport in the ab- 
sence of the reaction term. Generalizations of Fisher's 
equation that take into account finite correlation time 
effects in the transport term have been studied recently 

The logistic term u(l — u) bears similarity to sinu be- 
tween its zeros. Replacement of the former term by the 
latter makes Eq. (|^) the sine-Gordon equation which de- 
scribes, for example, the continuous limit of a chain of 
harmonically coupled pendulums The main differ- 
ence between the two equations is the periodicity of sin u 
and, consequently, an infinite number of zeros (equilib- 
ria). The logistic term provides only two equilibria. Ob- 
viously, we expect low amplitude solutions of Eq. (||) to 
be qualitatively similar to sine-Gordon's solutions, and 
large amplitude solutions (in particular solitary wave so- 
lutions connecting the two equilibria) to be quite different 
from the solitonic solutions of the sine-Gordon equation. 

We analyze in this paper traveling wave solutions of 
Eqs. ([|), which can be calculated exactly by analytical 
means. It is customary to analyze a piecewise linear sim- 
plification of a nonlinear system like Eq. (|3|), such as is 
done in Ref. Our analytic solution will provide us 
an opportunity to assess the validity of a piecewise lin- 
earization of the problem. Traveling waves of a related 
equation are analyzed in Section III. 



II. TRAVELING WAVES 
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We look for traveling- wave solutions of Eq. (g), moving 
in the direction of increasing x, by means of the ansatz: 
u(x — ct) = U{z). In general c, the speed of the nonlinear 
waves, will be different from v, the speed of the linear 



2 



waves dictated by the medium. The ordinary differential 
equation describing the shape of the wave is: 



(c 2 -v 2 )U" = kf(U) = -^. 



(4) 



A mechanical analogy involving a particle of mass m = 
|c 2 — v 2 \ is immediately available (see e.g. 0). Explicit 
solutions can be found by multiplying Eq. (0) by U' and 
integrating to find: 

l -(c 2 - v 2 )U' 2 -E = J kf(U)dU = -V(U), (5) 

where E is the "energy" of the particle, that depends on 
the initial conditions. Integrating once more: 



20 



dU' 



u ^-^[E-V(W)} 



(6) 



This is an elliptic integral that can be solved exactly 
since the radicand is a cubic function of U. This is done 



in Section II A 
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FIG. 1: The potentials of the mechanical analog of Eq. (|^). 
The cases shown are V and V , corresponding to faster than 
v and slower than v traveling waves respectively. 



this value. Since we are interested in situations where 
u describes a density or a population, we will exclude 
these solutions that contain a negative amplitude in the 
present discussion. Nevertheless, they can be analyzed 
without difficulty in the same way. 



A. General solutions 

As we noted, U(z) can be found exactly by perform- 
ing the elliptic integral in Eq. (H). Note the following 
simplified form H: 



dU 



VI - U 3 + U 2 



= AT arg sin VB + CU\j 



(7) 



where T((j)\^f) is the elliptic integral of the first kind, 
and A, B, C, and 7 are constants involving the roots 
of 1 — U 3 + U 2 . The constant 7 is termed the modulus of 
the elliptic integral, not to be confused with the parame- 
ter which is m = 7 2 . (We are following the nomenclature 
of Ref. ||, chapter 8.1, but denoting the modulus 7 in- 
stead of k to avoid confusion with the reaction constant.) 
The relation between the elliptic integral T and the jaco- 
bian elliptic function sn(z\j) suggests that the solution 
of Eq. (H) is of the form: 



U{z) = A + Bsn 2 (/3^|7), 



(8) 



where A, B, (3 and 7 are constants, depending on the ini- 
tial conditions, that remain to be determined. Differenti- 
ating Eq. ^ with respect to z, squaring and multiplying 
by m/2, we arrive at: 



-U' 2 = -2mpB(U-A) 



1 - W L (U-A) + ^L(U-A) 2 



(9) 



Equation (^|) is the kinetic energy that we can equate to 
E — V(U), which in turn we can write in the form: 



E - V{U) 



E 



U 2 



The potential of this mechanical system is V(U) = 
k(U 3 /3 ~ U 2 /2) if c > v, and V(U) = -V(U) if c < v 
(refer to Fig. [I] for a comparison) . The potential V has 
a maximum V = at U = and a minimum at U =1. 
Thus, if c > v, the system performs oscillations around 
[7 = 1, with the amplitude U always remaining positive, 
provided that the energy satisfies — k/6 < E < 0. If 
E = there is a solitary-wave solution with the shape of 
a traveling pulse, growing from U = to U = 1.5 and 
back to U = 0. 

The potential V, instead, has a minimum at U = and 
a maximum at U = 1 . The system oscillates around U = 
if the energy is positive and smaller than a finite value 
(0 < E < k/Q), and a traveling trough if it is exactly 



= -(U-Ui){U-Ui)(U-U z ), (10) 

where U\, U2 and U3 are the zeroes of E — V, and de- 
pend on the initial conditions. By equating Eq. (|l^) and 
Eq. (|^) we can express the unknowns A, B, (3 and 7 as 
functions of U\, U2, U3 and k. Because these expressions 
are intricate and do not add much to the analysis, we 
do not display them. However, one of the consequences 
of these expressions is the following dispersion relation, 
relating the speed of the nonlinear waves to system pa- 
rameters: 



^(a+V^-i^-Uz) 2 



-1 1/2 



(11) 
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FIG. 2: Half waves of the periodic solutions of different en- 
ergy. Greater amplitude (and lower frequency) correspond to 
higher energy. The parameters are k — 1, m = 2. 



where a = —U\ + 2[7 2 — U3. 

Equation (Q) can also be integrated numerically to 
yield z(U) for a given value of the energy E. In Fig. || 
we show half- waves for several energies. The complete 
solutions are periodic functions of z, symmetric around 
their maxima, and we are plotting here only one half of 
a period. These solutions correspond to waves of pro- 
gressively lower energy, smaller amplitude and, as can be 
seen, corresponding higher frequency. A grey line relat- 
ing amplitude to period connects the half periods, and 
can be seen to converge to a value of the half period of 
s/2tt ps 4.44, corresponding to the harmonic small oscil- 
lations. Waves of diverging period eventually become a 
traveling pulse, that we describe below. 



B. Validity of piecewise linearization procedures 

A piecewise linearization of Eq. (^) can be made, as 
in related models of reaction-diffusion processes (see for 
example @, ||, || for excitable systems, for an elec- 
trothermal instability, and the generalization of Fisher's 
equation found in ||). Following ||, the potential of the 
piecewise linear oscillator can be written as: 



-V kU 2 /{2a), 



U <a 



V ^ = <v"(UV2- bU ), U>a. < 12 > 



with free parameters, a, b and Vq, that can be adjusted 
to match the nonlinear one. The corresponding oscillator 
system is: 



mU" = —U (U<a), 



(13) 



whose solution can be found explicitly. Equations jl^ ) 
and ( |l4| ) describe harmonic oscillators, the first one with 
a positive restitutive force, so that the solution will be 
unstable (if m > 0). If the initial condition Uo is lower 
than a and has zero velocity, the general solution is a 
smooth match (a continuous match of U and U') of the 
solutions of both e qua tions ([l3 14). If Uq > a the solu- 



qua t 

s 



tion is just that of 

We define natural frequencies: u>i = y/k/m and u>2 = 
y/k/[m(b-a)]. The solution to (HJl^) is: 



U(z) 



Uq cosh(wiz), 

Acos(uj2z) + Bs'm(uj2z) + b, 



V>1 < 15 > 



where A and B are found as a solution to the following 
linear system: 

cos(w2Z a ) A + sin(uj2Z a ) B = a — b, 
-uj 2 sm(uj2Z a ) A + lo 2 cos(uj 2 z a ) B = Uqlux sinh(wxz a ), 

(16) 

where z a is the time it takes the unstable solution to grow 
from Uq to a, namely: z a = lu^ 1 argcosh(a/J7o). 

When the initial condition is greater than a, the solu- 
tion is: 



U{z) = {Uq - b) cos{uj2z) + b, (U > a). 



(17) 



mU" 



m(b — a) 



(U-b) (U>a), (14) 



The two kinds of solutions (those of the nonlinear sys- 
tem and those of the piecewise linear) have the same qual- 
itative shape. The piecewise linearization of the system 
is harmonic around the minimum at U — b, so there are 
harmonic oscillations of amplitude up to b— a (oscillations 
of greater amplitude feel the nonlinearity at U = a). The 
frequency of these oscillations is u>u n — yjk/[m(b — a)]. 
If 6=1, a = 1/2 and Vq = 1/4 (to match the logistic 
case U(l — U) as well as possible) then u>u n — Wk/2m 
which is \[2 smaller than what is found for the nonlinear 
system. 

This discrepancy can be seen in Fig. ^ where the re- 
lationship between the amplitude and the half period is 
shown, in a way similar to a dispersion relation. The 
nonlinear solution is shown as the full line. The piece- 
wise linear solution that preserves the intensity of the 
nonlinearity is shown as the dashed line. It's frequency 
at small amplitude overestimates the nonlinear one by a 
factor of v2- 

Of course, the piecewise linear procedure can be ad- 
justed to fit the parameter we want, such as the slope at 
[7 = 1. In this case the frequency of oscillations is \fkjm, 
but f(a = 1/2) = 1/2 instead of 1/4. This case is shown 
in Fig. H as the dotted line; it has the correct frequency 
at small amplitude. Both linearizations overestimate the 
amplitude of the nonlinear waves. Fig. |J shows the three 
potentials. In can be seen that the amplitude of traveling 
waves of a certain energy will always be overestimated in 
the piecewise linear system. 



4 



1.6 



1.2 - 



o 

3 

-*-» 

^ 0.: 

On 



0.4 



0.0 




purely harmonic regime 



nonlinear 
piecewise linear 
piecewise linear 
(preserving slope at t/=l) 
_i i i i i_ 



10 



12 



772 



FIG. 3: Amplitude vs. half period of traveling wave solutions. 
Results are shown for the fully nonlinear system and for two 
piecewise linearizations: one that preserves the value f(a) and 
one that preserves the slope of / at u = b, as shown in the 
legend. The parameters are: k = 1, m = 2, a = 1/2, 6 = 1. 



The speed of the pulse is related to its width. From the 
second moment of ( |l8| ) we find the following expression 
for the velocity as a function of the parameters: 



c(a) 



k ( V2tt 2 <j 2 



2/3 



(19) 



where a is the variance of ( |18| ) . It can be seen that wider 
pulses move faster than narrow ones, much in the way 
that the less steep fronts are faster in Fisher's equation 

0- 

This solitary wave solution is different from the kink- 
type soliton of the sine-Gordon equation, despite the sim- 
ilarity between the two equations. The potential in the 
sine-Gordon case is sinusoidal, allowing the state of the 
field u to shift by amounts of 2-7T in a kink. In the present 
case the potential V(u) continues growing as u — > oo, and 
the state u is forced back to u = after the excursion 
represented by the pulse (|I§|). 



III. FRONT SOLUTIONS FOR A RELATED 
SYSTEM 
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FIG. 4: The potential of the mechanical analog of the fully 
nonlinear and the piecewise linear systems. The case shown 
corresponds to c > v. 



C. Traveling pulse at E — 

When E = (and c > v) the solution is a traveling 
solitary pulse instead of a periodic wave. It corresponds 
to the case 7 = 1, when sn(</>|l) = tanh(0), and it repre- 
sents the trajectory of the particle leaving the maximum 
of V at U — towards the right (see Fig. ||) and bounc- 
ing back to U — after reaching the maximum amplitude 
U = 1.5. In this case, the integral (||) simplifies consid- 
erably, and the pulse can be written down explicitly as: 



U(z) 



tanh 



(18) 



It is possible to approach our model of coupled reac- 
tors from another direction. Suppose that at each site of 
the lattice we have a system described by a first order 
equation of chemical or population dynamics: 



du n 



(20) 



So, if / is logistic, the "force" (deriving (|20j) with respect 
to time) is the cubic f'{u)f(u) — k(\ — 2u)u{\ — u) and 
the corresponding generalization of the wave equation is: 

Qt2 = v 2 ^ + *«(l-u)(l-2«). (21) 

Equation (Elj) is very similar to the 4 equation of parti- 
cle physics M. But in Eq. ( |2l| ) the u A potential is inverted 
and it has as a consequence only one stable vacuum state 
instead of two. 

Solutions of Eq. (|2l]) are qualitatively similar to those 
of Eq. (|^) around and near the equilibrium u = 0. Solu- 
tions far from the equilibrium can be analyzed with the 
same methods as those of Eq. (^|) , since the integrals in- 
volved in (JgJ) are also elliptic. In fact, Eq. (^lj) is similar 
to a Klein-Gordon c/> 4 model, but with an inverted poten- 
tial. The solitary wave in this case is of the kink type, 
and it is straightforward to calculate its shape as z(U): 



z(U) 



m U(l-U)lag[(l-U)/U] 



v / U±/2-U 3 + U 2 /2 



(22) 



This solution is a traveling front, connecting the hyper- 
bolic points located at the two maxima of the potential, 
at U = and U = 1. It represents the invasion of one of 
the states by the other. A kink from U = to U = 1 as 
well as an "anti-kink" from U = 1 to U = are possible, 
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FIG. 5: Half waves of the periodic solutions of Eq. ( |2l[ ) (thin 
black lines) and kink solution (^) (thick grey line). Different 
lines correspond to different negative energies of the mechan- 
ical analog. The kink is the solution with zero energy. 

opening interesting possibilities of multiple interacting 
kinks, which will be analyzed elsewhere. 

Oscillatory wave solutions are plotted in Fig. [| for sev- 
eral initial conditions. Lower energies tend to accumu- 
late to the kink solution (p2|), shown with a thicker line 
in grey. 

IV. CONCLUSION 

We have analyzed two systems of reacting elements 
harmonically coupled to first neighbors. The equation of 



motion in the continuous limit is, in each case, a wave 
equation with a nonlinear reaction term. The transport 
character in the systems is, thus, fully coherent, in con- 
trast to a case such as that of Fisher's equation jjj, in 
which it is fully incoherent. Intermediate transport co- 
herence has been analyzed recently [|[ ||] through the in- 
corporation of memory functions and it has been shown 
how the incoherent (Fisher) limit may be obtained for 
a memory with infinitely fast decay. The opposite limit 
wherein the memory is constant is Eq. (|^) in the present 
paper, while a related system is Eq. (|2l|). The latter rep- 
resents a set of logistic systems, coupled harmonically. 
There is an interesting difference between its solutions 
and those of other extended reaction models. In Fisher's 
reaction-diffusion equation, traveling front solutions of a 
single kind are found. In the generalized system stud- 
ied by us elsewhere || , two kinds of fronts — a front and 
an "anti- front" — have been found, which however cannot 
coexist: each exists for a different set of system parame- 
ters. By contrast, in the purely wave-like system we have 
studied in the present paper, the coexistence of kinks and 
anti-kinks is certainly possible. This opens an interesting 
possibility of multiple interacting fronts. 
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